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We calculated, within the NRQCD factorization formalism, the leading color-octet con- 
tributions to if) production in photon-nucleon and electron-nucleon collisions. The ex- 
pressions obtained depend on the NRQCD matrix elements (Og ( 1 So)) and (Og( 3 Pj)). 
These matrix elements can be determined by fitting to experimental data. The color- 
octet contribution to if) photoproduction is in the forward region of phase space, where 
there may be large corrections to the NRQCD result from higher twist terms. As to 
ip leptoproduction we point out that the theoretical uncertainties plaguing the photo- 
production calculation vanish in the large momentum transfer limit. In this region of 
phase space the NRQCD formalism should be valid, making leptoproduction an ideal 
laboratory for testing the theory. 

1. Introduction 

The nonrelativistic QCD (NRQCD) factorization formalism developed by Bod- 
win, Braaten, and Lepage provides a rigorous theoretical framework within which 
quarkonium production can be studied El. A central result of this formalism is that 
the cross section for the inclusive production of a quarkonium state H is a sum of 
products having the form 

a(A + B - H + X) JpL^Q*) , (1) 

n m Q 

where vhq is the mass of the heavy quark Q 7 F n are short-distance coefficients, and 
are NRQCD four-fermion production operators with naive energy dimension 
d n . The short-distance coefficients, F ni are associated with the production of a 
QQ pair with quantum numbers indexed by n (angular momentum 2S+1 Lj and 
color 1 or 8). They can be calculated using perturbative techniques. The NRQCD 
production matrix elements, (O^) = (0|C^|0), parameterize the hadronization into 
H of a QQ pair with quantum numbers indexed by n. They can be determined 
phenomenologically. 



*Presented at the 1996 University of Illinois at Chicago Quarkonium Physics Workshop. This work 
was done in collaboration with Jim Amundson, Ivan Maksymyk, and Tom Mehen. 
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The power of the NRQCD formalism stems from the fact that Eq. (|l|) is es- 
sentially an expansion in the small parameter v 2 , where v is the average relative 
velocity of the O and Q in the boundstate H . v 2 ~ 0.3 for charmonium. NRQCD 
w-scaling rules □ allow one to estimate the relative sizes of the various (O^)- This in- 
formation, along with knowlege of the dependence of the F n on coupling constants, 
permits one to decide which terms must be retained in expressions for observables 
to reach a a given level of accuracy. Generally, to leading order, factorization for- 
mulas involve only a few matrix elements, so several observables can be related by 
a small number of parameters. Thus it is possible to test the NRQCD factorization 
formalism by determining if a body of data can be consistently fit by the most 
important NRQCD production matrix elements. Moreover, by studying the regime 
in which the NRQCD factorization formalism fails we can gain valuable insight into 
the limitations of the theory. 

As to J/ip production, in many instances the most important NRQCD matrix 
elements are (Of( 3 5i)), <0^( 8 Si)>, (Of^So)}, and (Of( 3 Pj)}. These matrix ele- 
ments appear in the expressions of rates for ip production in hadronic collisions, in 
Z° decayJn e + e~ annihilation, in photon-nucleon collisions, and in lepton-nucleon 
collisions u. Fitting the leading order predictions to the various experimental data 
has revealed some inconsistencies □ U. In particular thevalue determined for the lin- 
ear combination (Og ( 1 5'o)} + 3(0g ( 3 Po))/m^ at CDF □ appears to be incompatible 
with the value determined for the linear combination {Of ( 1 5o)) + 7 (Of ( 3 Po)) / m 2 
from photoproduction □ and other hadroproduction El experiments. Another as- 
pect of this problem is that if one uses the CDF measurement to make an estimate of 
the mag nitude of (Of^So)) and (Ot( 3 Pj)), assuming both matrix elements to 
positive, the color-octet contribution to inelastic tp photoproduction is too large I 
Throughout the paper we will refer to this inconsistency as the "photoproduction 
conundrum" . 

In this work we focus on aspects of ip photoproduction and electroproduction. 
Within the context of these processes we investigate some of the issues involved in 
applying the NRQCD factorization formalism, and study the possible limitations of 
the theory. More precisely our goal is to see if an analysis of ip leptoproduction can 
resolve the photoproduction conundrum. 

The paper is divided into two parts. The first part, section 2, is a calculation of 
the forward ip photoproduction rate. We fit the theoretical expression for the rate 
to experimental data and determine a value for the linear combination (Og ( 1 S'o)) + 
7(Og( 3 Pj))/m 2 . We point out that there are large corrections to this result; a sign 
that this process is not amenable to analysis within the NRQCD formalism. 

In the second part of the paper, section 3, we present a calculation of the ip 
leptoproduction rate. For large momentum-transfer squared, Q 2 , this process does 
not suffer the sizable corrections afflicting the photoproduction calculation, and is, 
therefore, well within the regime of applicability of the NRQCD formalism. Unfor- 
tunately, at this time, there is no data available on forward ip leptoproduction at 
large Q 2 . Thus we fit the theoretical result to experimental data covering low to 
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moderate values of Q 2 , and measure (OtfSo)) and (Ot( 3 P )). We find a negative 
value for (Of( 3 P )). By studying the rcnormalization of the P-wave operator we 
argue that this is not shocking. The values of the color-octet matrix elements deter- 
mined here resolve one aspect of the photoproduction conundrum: the CDF analysis 
is now consistent with photoproduction and other hadroproduction analyses. It is 
unclear, however, whether the other aspect of the photoproduction conundrum is 
resolved: that the color-octet contribution to inelastic ip photoproduction is too 
large. 

2. Photoproduction of ip 

One of the earliest calculations of ip production was carried out by Berger and 
Jones in 1981 In this paper the authors present a calculation of the rate for ip 
production in 7-nucleon collisions carried out within the color-singlet model. In 
this model one assumes that there is a nonzero probability for a QQ pair to form a 
hadron H only if the QQ pair is produced at short-distances (i.e. at distances on 
theorder of 1/ wiq or less) with the quantum numbers of the dominant Fock state of 
H E3. For example a cc pair has a nonzero probability of forming a ip only if the cc 
pair is produced at short-distances in a color-singlet state with angular-momentum 
quantum numbers Si. 

In order to gain a deeper understanding of the color-singlet results we need 
to define the parameter z = N ■ P/N ■ k. Here N is the initial-state nucleon 
four-momentum, k is the initial-state photon four-momentum, and P is the ip four- 
momentum. In the rest frame of the nucleon z is the fraction of the photon energy 
that is carried away by the ip. The region of phase space where z < 0.9 is, by 
convention, defined as "inelastic". The remaining region of phase space 0.9 < z < 
1.0 is defined as "forward". 

If one includes next-to-leading order QCD corrections and next-to-leading order 
rclativistic corrections the color-singlet model explains the experimental data in the 
inelastic region til. However, the color-singlet modeLprediction falls nearly an order 
of magnitude below the data in the forward regions. This is not unexpected since, 
as Berger and Jones pointed out over a decade ago in the forward region of phase 
space the final state gluon couples to on-shcll quark lines; thus a s for the vertex is 
ill-defined and the emission of the gluon is nonperturbative. 

The failure of the color-singlet model to explain forward ip photoproduction 
leads us to the question: can ip photoproduction be understood in the NRQCD 
factorization formalism? 

According to the NRQCD factorization formalism the photoproduction cross 
section is given by Eq. [I| The w-scaling rules tell us that the probability for produc- 
ing H from a QQ pair that does not have the quantum numbers of the dominant 
Fock state of H is down by powers of v 2 relative to the probability for producing 
H from a QQ pair with the quantum numbers of the dominant Fock state of H . 
As pointed out previously, v 2 is not in general zero so it is possible for a QQ pair 
produced at short-distances with any quantum numbers to hadronize into H . Only 
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Fig. 1. Leading order diagrams for the photoproduction of a cc in a color-singlet, 3 Si state. 




Fig. 2. Leading order diagrams for the photoproduction of a cc in a color-octet state with angular 
momentum configuration 1 So, 3 Poi an d 3 P2- 

in the v — ► limit is the color-singlet model recovered. 

Within the NRQCD factorization formalism the leading contributions to the 
photoproduction cross section come from the production of a cc pair in a color- 
singlet 3 S\ state, or from the production of a cc pair in a color-octet state in 
either a Sq or 3 Pj configuration. To leading order in a s the color-singlet short- 
distance coefficient Fi( 3 S'i) can be determined from the diagrams in figure 1. The 
result is proportional to the expression derived for the ip photoproduction cross 
section in the color-singlet model To leading order in a s the color-octet short- 
distance coefficients Fg( 1 5o) and F$( 3 Pj) can be determined from the Feynman 
diagrams in figure 2. The color-octet matrix elements are suppressed by v A w 
0.1 relative to the color-singlet matrix element, but the color-octet short-distance 
coefficients are enhanced by a n / a s (2m c ) f» 10 relative to the color-singlet short- 
distance coefficient. Thus both contributions are equally important. 

The NRQCD factorization formalism separates effects of short-distance scales of 
order m c or higher, which are associated with the production of the cc, from long- 
distance scales such as m c v, which are associated with the hadronization of the cc 
into the ip. This separation is embodied in the factored form of the production cross 
section as shown in Eq. ([[]). The spirit of the formalism is that 7 + g — > ccii 3 Si) + g 
is the correct short-distance process as long as the emission of the final state gluon 
takes place within a distance of order l/m c of the interaction point. This means 
that this gluon must have momentum of order m c or greater. If the final state gluon 
has momentum less than m c its emission is not part of the short-distance process 
and the diagrams shown in figure 1 do not describe the hard scattering. 

Therefore, we hypothesize that the color-singlet calculation is not valid for the 
region of phase space 0.9 < z < 1.0, since in this region the emission of the final 
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state gluon is nonperturbative, and is thus not part of the short distance process. 
Rather photoproduction in the forward region is described by the leading color- 
octet process. Inspired by the NRQCD formalism we will limit the color-singlet 
contribution to the region of phase space where < z < 1 — A, where A is some 
arbitrary cutoff of order v 2 . Then the color-octet contribution produces ip m the 
region 1 — A < z < 1. By convention in experiments the cutoff point between 
inelastic and forward ip production is chosen to be z w 0.9. Thus we choose A w 0.1. 
We wish to emphasize the cutoff does not arise naturally in the NRQCD formalism, 
rather it is an assumption made in the spirit of the formalism. 

The color-singlet contribution to ip photoproduction has been studied exten- 
sively Ed. Let us now consider the color-octet contribution. 

The leading color-octet contribution to the ip photoproduction cross section can 
be_calculated from the Feynman diagrams given in figure 2. The resulting expression 



CT ( 7 + N _> 4 + X) = fdx f g/N (x) as(2TOc) " e ^ 3 5(xs - Am 2 c ) © . (2) 
J m£ 

where x is the momentum fraction of the incoming gluon relative to the nucleon, 
fg/if(x) is the gluon distribution function for the nucleon, and 

G = {OteSo)) + 7 { °* (3 2 Po)) . (3) 
m c 

Since the values of the color-octet matrix elements are not known we can not make 
a prediction for the forward photoproduction cross section. However, we can fit the 
results to experimental data and make a prediction for ip production in some other 
process. 

Figure 3 shows a fit of Eq. (]2 ) to forward cross^section measurements from the 
fixed target experiments E687 B, NA14 M, E401 0, NMC El, and E516 0. Using 
a s (2m c ) = 0.26 and m c = 1.5 GeV we obtain the value 

9 = 0.02 GeV 3 , (4) 



No theoretical error has been quoted here. The expression given in Eq. (^|) is very 
sensitive to the value we choose for the parameter m c , which results in a large 
uncertainty in the value determined for 0. Given this error, the number presented 
in Eq. (^) is consistent with the value = 0.03 GeV 3 measured in nN collisions B. 

There is an important point to make regarding the calculation of forward ip 
photoproduction. Namely, in the derivation of the factorization formula presented 
in Eq (Q) higher twist terms have been neglected. In the forward region these higher 
twist terms can be large E3. For example, they could describe correlations between 
the initial and final state resulting in diffractive and elastic processes. Such processes 
can have large contributions to the ip forward photoproduction cross section, and 
can, therefore, significantly affect the value determined for Ea. 
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Fig. 3. A one parameter fit to experimental data of the NRQCD factorization formalism result for 
forward photoproduction of ip. 



Given the large theoretical uncertainty associated with the forward photopro- 
duction calculation it is best to regard the value presented in Eq. as an upper 
limit. In fact, if the value determined for in photoproduction is consistent with 
fits of the color-octet matrix elements to other ip production data we can be certain 
that higher twist corrections, and diffractive and elastic contributions are small. 

3. Leptoproduction of ip 

Many of the theoretical uncertainties plaguing the photoproduction calculation can 
be avoided by requiring the incoming photon to be highly virtual. This introduces 
a new scale into the process: Q 2 the momentum transfcrcd through the photon 
squared. Higher twist terms will be suppressed by powers of Q 2 , and will, therefore, 
vanish in the large Q 2 limit. 

A process in which the photon can be highly virtual is ip leptoproduction. The 
calculation of ip leptoproduction is analogous to the photoproduction calculation, 
except the incoming photon is off shell. In the region of phase space where 0.9 < 
z* < 1.0 (the * denotes that the photon is virtual) the leading contribution comes 
from the fusion of the virtual photon and the gluon into a color-octet cc pair in 
either a 1 Sq, 3 Pq, or 3 P2 configuration, followed by the hadronization of the cc pair 
into a ip. The Feynman diagrams used to determine the short-distance coefficients 
are shown in figure 4. 

?he expression for the cross section determined from the diagrams in figure 4 



is 



s 

~Q 2 ~ J ~y J X f g/N ^ 5(yXys ~ 4m ^ 
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Fig. 4. A two parameter fit to EMC data of the NRQCD factorization formalism result for forward 
leptoproduction of ip. 



2a 8 (^)a 2 eiTr 2 1 + (1 - y)< 



xsm c y 
8(2m c ) 2 Q 2 (Ot( 3 Pj)) 



(otes )) 



3Q 2 + 7(2m c ) 2 (Ot( 3 Pj)) 



xys 



x 2 ys 2 



D)7. 



(5) 



where /x 2 = Q 2 + m 2 , and x and f g /N( x ) are the same as in Eq. (Q). The momentum 
fraction of the virtual photon relative to the incoming lepton is y ~ N ■ q/N ■ k, 
where N is the nucleon four-momentum, q is the photon four momentum, and k is 
the incoming lepton four-momentum. 

The result presented in Eq (||) holds for all values of Q 2 . Taking the limit Q 2 — > 
one recovers the photoproduction result convoluted with the electron splitting 
function: 



lim cr(e + P -> e + ib + X) 
Q 2 ^o 



a 
27 



dQ 2 

Q 2 



1 l + (l- y )a 
2/ 



As discussed previously, in this limit corrections to the cross section from higher 
twist terms, may be large. However, in the high-energy limit Q 2 ,s 3> (2m c ) 2 we 
expect contributions from higher twist terms to vanish. Letting Q 2 , s 3> (2m c ) 2 we 
obtain 



l/Q 



dy 



l + (l-y) 2 



lim a(e + P^e + yj + X)^^- [ [ 

2 ,m2/s->0 2tt J Q 2 J 

, , 4a s (Q 2 )ae 2 y ( ' ,^ flo ^ , {Oj( s Pj)) \ £ , 

dx f g/N (x) — I (Og( So)) + 3 — 2_ I s(xys 



J )-(7) 



Note that this expression does not depend very strongly on the value chosen for 
m c . Since theoretical corrections to Eq. (Q) are expected to be small, high energy 
leptoproduction provides an excellent means to measure the linear combination 
(Of? So)) + 3(C^( 3 Pj))/m 2 . This is precisely the linear combination of NRQCD 
matrix elements determined from CDF data on ip production at high transverse 
momentum. Therefore, ip leptoproduction data taken in the high energy limit will 
provide us with the opportunity to test the NRQCD factorization formalism by 
measuring the same linear combination of matrix elements in two different processes. 

As of yet there is no ifj leptoproduction data that truely falls in the high energy 
regime. Therefore we will naively fit Eq. (S) to the available data. The result of 
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Fig. 5. A two parameter fit to EMC data of the NRQCD factorization formalism result for 
forward leptoproduction of ip. The horizontal axis is Q 2 , and the vertical axis is da(-y* + P — » 
i> + X)/dlogQ 2 . 

a fit to EMC data over the entire range of Q 2 (0 < Q 2 < 15 Gev 2 ) is shown in 
figure 5. The values of the color-octet matrix elements determined from this fit are: 

(OtCS )) = 0.04 GeV 3 , 

(0 3 (3 f j)) = -0.003 GeV 3 . (8) 

Note that (Of( 3 Pj)) is negative! Does this mean that the NRQCD factorization 
formalism fails to describe tp leptoproduction? Perhaps; remember there could be 
large corrections to the low Q 2 region. However, there is another explanation. 

Loop corrections to NRQCD operators give rise to ultraviolet power divergences 
which have the form of renormalization of lower-dimensional operators. Since di- 
mension six is the lowest possible energy dimension of NRQCD production operators 
the dimension-six operator Of ( 1 S'o) will not have a power-divergent contribution. 
However, the dimension-eight operator 0g ( 3 Pj), will have power divergences pro- 
portional to dimension six operators. At order a Sl using a momentum cutoff A, 
these divergences can be removed by defining the renormalized operator to be 

Ot( 3 Pj)r = Ot( 3 P.j)o ~ c ias A 2 Of ( 3 S X ) - c & a s k 2 Ot( 3 S,) , (9) 

where the coefficients c\ and c% are adjusted to cancel quadratic divergences at 
order a s . Matrix elements of the bare operator Og ( 3 Pj)o must indeed be positive 
definite, but this is not necessarily true for the renormalized operator. If one defines 
the operator 0%{ 3 Pj) using a method like dimensional regularization which auto- 
matically removes power divergences, the regularization method makes the above 
subtractions implicitly. Therefore the renormalized operator need not have positive 
matrix elements til. 
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In particular it is easy to see that the most important subtracted term in Eq. (|9]) 
is the one proportional to the matrix element Of ( 3 Si). This term is suppressed by 
a factor of a s , but it is enhanced by a quadratic divergence and a relative factor of 
1/f 4 . Thus it is not too shocking if the resulting renormalized matrix element turns 
out to be negative. Note that renormalization respects the w-scaling rules which 
require the magnitude of the matrix element of 0g ( 3 Pj) to scale as v 7 . 

Accepting the explanation proffered above let us compare the results presented in 
Eq. (||) to photoproduction and hadroproduction results. The numbers determined 
in leptoproduction are consistent with the photoproduction result given in Eq. (Q), 
and are, therefore, consistent with the result of an analysis of TtN collisions E. Using 
CDF data on ip production at the Tevatron Cho and Leibovich have determined H 

V/Y3 i 



(OfCSo)) + 3 <0s( = 0.066 GeV 3 



(10) 



Substituting the values given in Eq. (g) into the left-hand side of Eq. (|1C|) we obtain 
0.03 GeV 3 . Given the large theoretical uncertainty associated with calculations of 
ip production in hadronic collisions these results are consistent. 

4. Conclusion 

We have calculated, within the NRQCD factorization formalism, the leading color- 
octet contributions to ip photoproduction and ip leptoproduction. The expressions 
obtained depend on the color-octet matrix elements (Of ( 1 5'o)) and (0g ( 3 P )). The 
NRQCD factorization formalism may be tested by fitting these matrix elements to 
experimental data, and then making predictions for ip production in some other 
process. 

The color-octet contribution to the photoproduction cross section is in the region 
of phase space where the ip is produced in the forward direction. In this region higher 
twist terms, which were neglected in the derivation of the factorization formula 
Eq. |l|, may be large. Therefore, forward ip photoproduction is not amenable to 
testing the NRQCD factorization formalism. However, we can learn about the 
limitations of the theory, and the size of possible corrections. 

The leading color-octet contribution to the ip leptoproduction cross section does 
not suffer from the same problems as the photoproduction calculation if we restrict 
ourselves to the high Q 2 regime. In this region ip leptoproduction should provide 
an ideal laboratory for testing the NRQCD factorization formalism. However, there 
is at this time no experimental data for ip production in the asymptotic regime. 
Therefore, keeping in mind that there may be large corrections to the cross section 
in the low Q region, we fit our results to EMC data for < Q 2 < 15 GeV 2 . We 
determine (O^^So)) = 0.04, and (O 8 A ( 3 P )} = -0.003. The negative value for 
the P-wave matrix element is acceptable since the renormalized P-wave operator 
is a subtraction of two divergent terms, Of ( 3 P/)r = Of ( 3 P/)o - c 1 a s A 2 Of ( 3 Si), 
where the second term on the right-hand-side can be larger than the first term 
on the right-hand-side. The values determined for the color-octet matrix elements 
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resolve one aspect of the photoproduction conundrum: the discrepancy between the 
CDF analysis and photoproduction and other hadroproduction analysis. However 
it is not clear if the other aspect of the photoproduction conundrum is resolved: the 
color-octet contribution to inelastic ip production is too large. 
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